Fractional topological liquids with time-reversal symmetry and their lattice realization 
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We present a class of time-reversal-symmetric fractional topological liquid states in two dimensions 
that support fractionalized excitations. These are incompressible liquids made of electrons, for which 
the charge Hall conductance vanishes and the spin Hall conductance needs not be quantized. We then 
analyze the stability of edge states in these two-dimensional topological fluids against localization 
by disorder. We find a Z 2 stability criterion for whether or not there exists a Kramers pair of edge 
modes that is robust against disorder. We also introduce an interacting electronic two-dimensional 
lattice model based on partially filled flattened bands of a Z 2 topological band insulator, which 
we study using numerical exact diagonalization. We show evidence for instances of the fractional 
topological liquid phase as well as for a time-reversal symmetry broken phase with a quantized 
(charge) Hall conductance in the phase diagram for this model. 



I. INTRODUCTION 

The hallmark of the integer quantum effect (IQHE) 
in an open geometry is the localized nature of all 
two-dimensional (bulk) states while an integer number 
of chiral edge states freely propagates along the one- 
dimensional boundaries"^. These chiral edge states are 
immune to the physics of Anderson localization as long 
as backward scattering between edge states of opposite 
chiralities is negligible^'. 

Many-body interactions among electrons can be 
treated perturbatively in the IQHE provided the charac- 
teristic many-body energy scale is less than the single- 
particle gap between Landau levels. This is not true 
anymore if the chemical potential lies within a Landau 
level as the non-interacting many-body ground state is 
then macroscopically degenerate. The lifting of this ex- 
tensive degeneracy by the many-body interactions is a 
non-perturbative effect. At some "magic" filling frac- 
tions that deliver the fractional quantum Hall effect 
(FQHEjP", a screened Coulomb interaction selects a 
finitely degenerate family of ground states, each of which 
describes a featureless liquid separated from excitations 
by an energy gap in a closed geometry. Such a ground 
state is called an incompressible fractional Hall liquid. 
The FQHE is an example of topological ordejP^. In 
an open geometry, there are branches of excitations that 
disperse across the spectral gap of the two-dimensional 
bulk, but these excitations are localized along the direc- 
tion normal to the boundary while they propagate freely 
along the boundarySSHUl Contrary to the IQHE, these 
excitations need not all share the same chirality. How- 
ever, they are nevertheless immune to the physics of An- 
derson localization provided scattering induced by the 
disorder between distinct edges in an open geometry is 
negligible. 

The integer quantum Hall effect (IQHE) is the 
archetype of a two-dimensional topological band insu- 
lator. The two-dimensional Z 2 topological band insula- 



tor is a close relative of the IQHE that occurs in semi- 
conductors with sufficiently large spin-orbit coupling but 
no breaking of time-reversal symmetryS^Hni As with the 
IQHE, the smoking gun for the Z 2 topological band insu- 
lator is the existence of gapless Kramers degenerate pairs 
of edge states that are delocalized along the boundaries of 
an open geometry as long as disorder-induced scattering 
between distinct boundaries is negligible. In contrast to 
the IQHE, it is the odd parity in the number of Kramers 
pairs of edge states that is robust to the physics of An- 
derson localization. 

A simple example of a two-dimensional Z 2 topological 
band insulator can be obtained by putting together two 
copies of an IQHE system with opposite chiralities for 
up and down spins. For instance, one could take two 
copies of Haldane's modeiPD each of which realizes an 
integer Hall effect on the honeycomb lattice, but with 
Hall conductance differing by a sign. In this case the spin 
current is conserved, a consequence of the independent 
conservation of the up and down currents, and the spin 
Hall conductance inherits its quantization from the IQHE 
of each spin species. This example thus realizes an integer 
quantum spin Hall effect (IQSHE). However, although 
simple, this example is not generic. The Z 2 topological 
band insulator does not necessarily have conserved spin 
currents, let alone quantized responses. 

Along the same line of reasoning, two copies of a FQHE 
system put together, again with opposite chiralities for 
up and down particles, would realize a fractional quan- 
tum spin Hall effect (FQSHE), as proposed by Bernevig 
and ZhangPSl (See also Refs.lT§landl2Tjl) Levin and Stern 
in Ref. 21 proposed to characterize two-dimensional frac- 
tional topological liquids supporting the FQSHE by the 
criterion that their edge states are stable against disorder 
provided that they do not break time-reversal symmetry 
spontaneously. 

In this paper, we shall not impose the condition that 
projection about some quantization axis of the electron 
spin from the underlying microscopic model is a good 



2 



quantum number. We will only demand that time- 
reversal symmetry holds. We shall thus distinguish the 
generic cases of fractional topological liquids with time- 
reversal symmetry from the special cases of fractional 
topological liquids with time-reversal symmetry and with 
residual spin-1/2 U(l) rotation symmetry. In the former 
cases, the electronic spin is not a good quantum number. 
In the latter cases, conservation of spin allows for the 
FQSHE. 

The subclass of incompressible time-reversal- 
symmetric liquids that we construct here is closely 
related to Abelian Chcrn-Simons theories. Other 
possibilities that are not discussed in this public ation, 
may include non-Abelian Chern Simons theorie^2122l 
or theories that include, additionally, conventional local 
order parameters (Higgs fields p2 

The relevant effective action for the Abelian Chern- 
Simons theory is of the forrrP^ 



S \ — Sq -\- S e -\- S s 



(1.1a) 



where 



S :=-J dtd 2 x e^> ^K tJ a* d v a? p , (1.1b) 



and 



S, 



S Q 



dtd 2 xe^—S t B^a p . 



(1.1c) 
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The indices i and j run from 1 to 2N and any pair thereof 
labels an integer-valued matrix element of the sym- 
metric and invertible 2N x 2N matrix K. The indices 
/i, is, and p run from to 2. They either label the compo- 
nent x M of the co-ordinates (t, x) in (2 + l)-dimensional 
space and time or the component A p (t,x) of an exter- 
nal electromagnetic gauge potential, or the component 
B^(t,x) of an external gauge potential that couples to 
the spin-1/2 degrees of freedom along some quantiza- 
tion axis, or the components of 2N flavors of dynamical 
Chern-Simons fields a^(i, x). The integer- valued compo- 
nent Q i of the 2iV-dimensional vector Q represents the 
i-th electric charge in units of the electronic charge e and 
obeys the compatibility condition 



R Q * = (-) 



K: 



(Lie) 



for any i = 1, ...,2N in order for bulk quasiparticlcs 
or, in an open geometry, quasiparticles on edges to obey 
a consistent statistics. The integer-valued component 
S i of the 2/V-dimensional vector S represents the i-th 
spin charge in units of the spin charge s along some 
conserved quantization axis. The operation of time re- 
versal maps A^(t,x) into A v (—t,x); B (t,x) into 
-g^ B v (-t,x); a^{t,x) into -g» u a l + N (-t, x) for i = 
1, . . . , TV and vice versa. Here, g = diag(+, — , — ) is the 



Lorentz metric in (2 + l)-dimensional space and time. We 
will show that time-reversal symmetry imposes that the 
matrix K is of the block form 



K = 



K 

A T 



A 

— K 



A T = —A , 



(1.2a) 



(1.2b) 



where k and A are N x N matrices, while the integer- 
charge vectors Q and S are of the block forms 



Q = 



s = 



Q 



(1.2c) 



The K matrix together with the charge vector Q and 
spin vector S that characterize the topological field the- 
ory with the action 1 1.1a I define the charge filling frac- 
tion, a rational number, 



Q T K" 1 Q 



(1.3a) 



and the spin filling fraction, another rational number 

1 



v a :=-Q' K-'S, 



(1.3b) 



respectively. The block forms of K and Q in Eq. (1.2) 
imply that 



= 0. 



(1.3c) 



The "zero charge filling fraction" ( 1.3c[ ) states nothing 
but the fact that there is no charge Hall conductance 
when time-reversal symmetry holds. On the other hand, 
time-reversal symmetry of the action (1.1a) is compatible 



with a non-vanishing FQSHE as measured by the non- 
vanishing quantized spin-Hall conductance 



2tt 



x 



(1.3d) 



The origin of the FQSHE in the action (1.1a 
U(l)xU(l) gauge symmetry when (2 



is the 
l)-dimensional 
space and time has the same topology as a manifold with- 
out boundary. We shall always assume that the U(l) 
symmetry associated with charge conservation holds in 
this paper. However, we shall not do the same with the 
U(l) symmetry responsible for the conservation of the 
"spin" quantum number. 

The special cases of the FQSHE treated in Refs. [T51 
and |^ correspond to imposing the condition 



A = 



(1.4) 



on the K matrix in Eq. (1.2a). This restriction is, how- 



ever, not necessary to treat either the FQSHE or the 
generic case when there is no residual spin-1/2 U(l) sym- 
metry in the underlying microscopic model. 

The effective topological field theor y (|1.1[ ) with the 
condition for time-reversal symmetry (1.2) is made of 
27V Abelian Chern-Simons fields. As is the case with 
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the FQHE, when two-dimensional space is a manifold of 
genus one without boundary (i.e., when two-dimensional 
space is topologically equivalent to a torus) , it is charac- 
terized by distinct topological sectors^^l. All topologi- 
cal sectors are in one-to-one correspondence with a finite 
number Af GS of topologically degenerate ground states 
of the underlying microscopic theorj0E3. This degener- 
acy is nothing but the magnitude of the determinant K in 
Eq. ( 1.1a I, which is, because of the block structure (1.2a|, 



in turn given by 



GS 



det 



K A 



/ A T -K 

\k A 

( A T — K 
U A. 
(integer) 2 . 



det 
Pf 



(1.5) 



To reach the last line we made used of the fact that the 
K matrix is integer valued. We thus predict that the 
class of two-dimensional time-reversal-symmetric frac- 
tional topological liquids, whose universal properties are 
captured by Eqs. (1.1) and (1.2), are characterized by 



a topological ground state degeneracy that is always the 
square of an integer, even if A ^ 0, when space is topolog- 
ically equivalent to a torus. (Notice that the condition 
that A is anti-symmetric implies that non-vanishing A 
can only occur for N > 1.) 

We discuss in detail the stability of the edge states as- 
sociated with the bulk Chern-Simons action ( 1.1 ) obey- 



ing the condition for the time-reversal symmetry (1.2) 



We consider a single one-dimensional edge and construct 
an interacting quantum field theory for 1 < N K < N 
pairs of Kramers degenerate electrons subject to strong 
disorder that preserves time-reversal symmetry. [The in- 
teger 2N K is the number of odd charges entering the 
charge vector Q (RefP^J)]. We identify the conditions 
under which at least one Kramers degenerate pair of elec- 
trons remains gapless in spite of the interactions and dis- 
order. Our approach is here inspired by the stability 
analysis of the edge states performed for the single-layer 
FQHE by Haldane in Ref. OH (see also Refs. [27] and [21]) , 
by Naud et al. in Refs. [29] and [30] for the bilayer FQHE, 
and specially that by Levin and Stern in Ref. f5T] for the 
FQSHE. As for the FQSHE, our analysis departs from the 
analysis of Haldane in that we impose time-reversal sym- 
metry. In this paper, we also depart from Ref. [2T]by con- 
sidering explicitly the effects of the off-diagonal elements 
A in the K matrix. Such terms are generically present 
for any realistic underlying microscopic model indepen- 
dently of whether this underlying microscopic model sup- 
ports or not the FQSHE. When considering the stability 
of the edge theory, we allow the residual spin- 1/2 U(l) 
symmetry responsible for the FQSHE to be broken by 
interactions among the edge modes or by a disorder po- 
tential. Hence, we seek a criterion for the stability of 
the edge theory that does not rely on the existence of a 



quantized spin Hall conductance in the bulk as was done 
in Ref. Ell 

The stability of the edge states against disorder hinges 
on whether the integer 



R := rg T (k - A)" 



(1.6) 



is odd (stable) or even (unstable). The vector g tog ethe r 
with the matrices k and A were defined in Eq. (1.2). 



The integer r is the smallest integer such that all the N 
components of the vector r (k — A) -1 g are integers. We 
can quickly check a few simple examples. First, observe 
that, in the limit A = 0, we recover the criterion derived 
in Ref. 1211 Second, when we impose a residual spin-1/2 
U(l) symmetry by appropriately restricting the interac- 
tions between edge channels, = —v^ = g T (k — A) -1 g 
can be interpreted as the Hall conductivity er xy in units 
of e 2 /h for each of the separately conserved spin com- 
ponents along the spin quantization axis. The integer r 
has the interpretation of the number of fluxes needed to 
pump a unit of charge, or the inverse of the "minimum 
charge" of Ref. 21 , Further restricting to the case when 
k = 1 N gives R = N (i.e., we have recovered the same 
criterion as for the two-dimensional non-interacting Z 2 
topological band insulator). 

When there is no residual spin-1/2 U(l) symmetry, 
one can no longer relate the index R to a physical spin 
Hall conductance. Nevertheless, the index R defined in 



Eq. ( 1.6 ) discriminates in all cases whether there is or not 



a remaining branch of gapless modes dispersing along the 
edge. 

Before we turn to the detailed analysis of the stability 
of the edge theory in Sec. Ill we shall first pose and 



answer the question of whether one can realize examples 
of the Abelian Chern-Simons subclass of time-reversal- 
symmetric topological spin states in a two-dimensional 
lattice model in Sec. [H] We construct extensions of the 
lattice models studied in Refs. STJ [32j [33] and [34] for 
which a FQHE was found by partially filling flat bands 
with non-trivial Chcrn numbers, as proposed in Refs. 1311 
l35| and [36l (See also Ref. [37] for a discussion of isolated 
flat bands with broken time-reversal symmetry in two- 
dimensional lattices; Refs. [551 and I3T)1 for recent progress 
on the understanding of the relations between Chern and 
Landau bands; and Ref. 40 for predicting that materials 
belonging to the family of heterostructures of transition- 
metal oxides, say LaAu0 3 , might realize time-reversal 
symmetric topologically non-trivial band insulators with 
nearly flat bands.) 

The systems studied here start with flat bands that 
realize at half-filling a two-dimensional integer quantum 
spin Hall band insulator. We study with the help of exact 
diagonalization the nature of the ground state selected by 
interactions at partial filling 2/3 of the lowest band. We 
find supporting evidences for a featureless ground state 
that is consistent with the existence of a spectral gap and 
a topological degeneracy 3 2 in the thermodynamic limit, 
associated with a A = state (i.e., a FQSHE driven by 
interactions in a region of the phase diagram) . This state 
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is unstable to spontaneous symmetry breaking of time- 
reversal symmetry induced by sufficiently strong interac- 
tions, which lands the system onto a state with degener- 
acy 3 (not the square of an integer), which we identify 
with a 1/3 FQHE of a magnetized state. 

We close this paper with a summary in Sec. |IV| We also 
include two appendices to render the paper reasonably 
self-contained. 



II. EXACT DIAGONALIZATION STUDY OF A 
TWO-DIMENSIONAL LATTICE MODEL WITH 
TIME-REVERSAL SYMMETRY 



In Sec 
question: 



III 



we will pose and answer the following 
Given a time-reversal-symmetric incompress- 
ible liquid-like ground state whose universal properties 
are encoded by a low-energy and long-wavelength effec- 
tive quantum field theory for 2N Abelian Chern-Simons 
fields, under what conditions is a Kramers pair of edge 
modes that propagates at the boundary protected against 
Anderson localization as long as time-reversal symmetry 
is preserved. In this section, we want to address the ques- 
tion, if and when such a posited topological state emerges 
in the first place. If an incompressible state is connected 
to a translation invariant band insulator, once interac- 
tions are switched off adiabatically, the answer is en- 
tirely governed by the Bloch states of the single-particle 
Hamiltonian and is well understoocP^^l. If, however, 
the incompressibility of the state emerges from the inter- 
actions, the problem is qualitatively different. 

As a starting point to study the second situation, we 
shall follow the approach of Refs. [15] and [21] and con- 
sider two decoupled copies of the same incompressible 
fractional quantum Hall (FQH) state to compose an in- 
compressible time-reversal symmetric fractional quantum 
spin Hall (FQSH) state from them. Then, we know that 
the ground state of the system is the direct product of 
the two FQH states and thus it is also an incompressible 
liquid. It is then natural to ask how stable this state is 
when the two FQH states are coupled, that is, whether in- 
teractions between the two FQH states are (i) destroying 
the incompressibility, (ii) breaking spontaneously time- 
reversal symmetry, or (iii) generating other incompress- 
ible time-reversal symmetric states that are not captured 
by the Abelian Chern-Simons theory in Eq. (1.1). 



We are going to address this question by numerical ex- 
act diagonalization of an interacting lattice model, where 
we will find evidence for all scenarios (i), (ii), and (iii). 
It is advantageous to consider a system in which the z 
component of the electronic spin is conserved [i.e., with 
a residual spin-1/2 U(l) symmetry]. This makes larger- 
system sizes accessible and allows to study the model 
with spin-dependent twisted boundary conditions to de- 
termine the ground state degeneracy as explained below. 
We consider spinfull electrons hopping on the square lat- 
tice A = A U B made up of the two sublattices A and B. 
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FIG. 1. (Color online) Numerical exact diagonalization re- 
sults of Hamiltonian ( |2.1[ ) for 16 electrons when sublattice 
A is made of 3 x 4 sites and with t 2 /t 1 — 0.4. (a) Ground 
state degeneracies. Denote with E n the n-th lowest energy 
eigenvalue of the many-body spectrum where E 1 is the many- 
body ground state (i.e., E n+1 > E n for n = 1,2,.. .). Define 
the parameter e n := (E n+1 — E n )/(E n — E ± ). If a large gap 
-EVi+i — E n opens up between two consecutive levels E n+1 and 
E n compared to the cumulative level splitting E n — E 1 be- 
tween the first n many-body eigenstates induced by finite-size 
effects, then the parameter e n is much larger than unity. The 
parameter e n has been evaluated for n — 3 and n = 9, yielding 
the red and blue regions, respectively. For all other n / 1, no 
regions with e n > O(l) of significant size were found. Within 
the limited range of available system sizes, it is thus not pos- 
sible to decide on whether and how the level-splitting above 
the ground state in the white regions of the parameter space 
extrapolates in the thermodynamic limit, (b)-(d) The lowest 
eigenvalues with spin-dependent twisted boundary conditions 
as a function of the twisting angle j x . The number of low- 
lying states that are energetically separated from the other 
states is 9, 3, and 3, respectively. In panel (c), it is the lowest 
band parametrized by j a that is three-fold degenerate. 



The Hamiltonian 



H = H n 



(2.1) 



decomposes into a quadratic part H a and an interacting 
contribution H int . 

First, let us define H , which consists of two copies 
of the 7r-flux phase with flat bands that was studied in 
Ref. [3ll one copy for each spin-1/2 species. We denote 
with cl the creation operator for an electron with 
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lattice momentum k and spin a =~f, ^ in the sublattice 
a = A, B and combine them in the sublattice-spinor 

ipl, a '■= ( c l a o-' c l b o-) ■ Then, the second quantized 
single-particle Hamiltonian reads 



H := 



B 



.,1,7 \ l^fcl 



B 



fceBZ 



Lb 




where the three-vector B k is defined by 
B a.k '■= °, 



+ tj e +i7r/4 (e" ife « + e 



3,fc 



2i 2 (cos k x — cos fc 



(2.2b) 

(2.2c) 
(2. 2d) 



and the three Pauli matrices r = (t x , r 2 , t 3 ) act on the 
sublattice index. Here, t 1 and i 2 represent the nearest 
neighbor (NN) and next-nearest neighbor (NNN) hop- 
ping amplitudes. The Hamiltonian (2.2a) is only well 
defined if t 2 ^ 0. 

One verifies that H is both time-reversal symmetric 
and invariant under spin-1/2 U(l) rotations. Indeed, the 
time-reversal operation T acts on numbers as complex 
conjugation and on the electron-operators as 



<4 



T 



fe,t 



fe,4.) 



4 



(2.3) 



The action of the spin-1/2 U(l) rotation TZ 1 by the angle 
< 7 < 2tt is given by 



fe.t' 



(2.4) 



The spectrum of H is gaped and comprises four dis- 
persionless bands with the energy eigenvalues 



£k,a.± — ±1 5 



cr=t4 



(2.5) 



Denoting the corresponding single-particle eigenstates by 
X fecr± , cr =tj4-) wc can define the spin-resolved first 
Chern number for each of the two pairs of degenerate 
single-particle bands as 



1 f d 2 /s , t _ 

C s,± ~ 2 / 2id Vk A Ufe,t,± V fe^fe,t,± 



fceBZ 



(2.6) 



We find G 



s.± 



±1. As a consequence, the non- 
interacting model exhibits an IQSHE with spin-Hall con- 
ductivity o~? = 2 x e/(47r) if the chemical potential lies 
in the single-particle spectral gap. 

The repulsive interactions in this model are defined by 



(2.7) 



where (zj) are directed NN bonds of the square lattice 
A = A U B, p ia is the occupation number of the site 
i with electrons of spin a. The interacting Hamilto- 
nian H lnt comprises an on-site Hubbard term with the 
coupling U > and a NN term which is parametrized 
by the coupling V > and the dimensionless number 
A € [0,1]. The value A = 1 corresponds to the spin- 
1/2 SU(2)-symmetric limit, while all other values of A 
correspond to the spin-1/2 U(l)-symmetric limit. These 
interactions lift the macroscopic degeneracy of the single- 
particle bands. They couple the spin- up and the spin- 
down sectors, if at least one of U or A is non-vanishing. 
Notice that H int shares both the time-reversal and spin- 
1/2 U(l) symmetries of the single-particle Hamiltonian 
H . 

Periodic boundary conditions are imposed on lattice 
A whereby sublattice A contains L x x L y sites. We fix 
the number of electrons to be 16 while L x = 3 and 
L = 4. We define the filling fraction 2/3 to be the 
number of particles, 16, divided by the number of Bloch 
single-particle states in the lowest spin-degenerate band, 
(3 x 4 x 2 x 2) /2 = 48/2 = 24. We then project Hamilto- 
nian (2.1 ) onto the states in the two lower single-particle 
bands Efe.f,-, Efe,j.,-i thereby assuming that the single- 
particle gap is much larger than the energy scale of the in- 
teractions. Exact diagonalization yields the many-body 
spectrum as a function of the interaction parameters 
A and U/V. We identify three distinct incompressible 
states in the X-U phase diagram [see Fig. [ija)]. 

Case X = U/V = 0: decoupled FQH states. The model 
decouples into two FQH-like states at 2/3 filling, one for 
each spin orientation. The low-energy effective theory for 
this state could be compatible with the choice 



K 



+1+10 

+1-2 
0-1-1 

0-1 +2/ 



Q 



(2. 



for the K matrix and the charge vector Q in that it has 
degeneracy |det-ftT| = 3 2 = 9 as confirmed by the nu- 
merical results. This phase is destabilized by introduc- 
ing a sufficiently strong coupling between the two FQH 
states via A and U . The ability to make a quantitative 
statement on the boundary of this phase in the phase 
diagram is here limited by the difficulty in deciding on 
the compressibility of a state from extrapolation of ex- 
act diagonalization studies of small systems sizes to the 
thermodynamic limit. 

Case X = 1, U/V > 2: Spontaneous symmetry break- 
ing. We observe that the ground state has the maximal 
spin-polarization that is allowed by the Pauli principle. 
To interpret this numerical result, first recall that, af- 
ter projection onto the lowest bands, at most L x x L y 
electrons may have the same spin (i.e., 12 for the case 
at hand). Now, the filling fraction is 2/3 (i.e., there are 
4/3xL x xL y = 16 electrons). If 12 electrons are fully spin 
polarized, which is what we observe numerically, then the 
remaining 1/3 x L x x L y = 4 electrons may form a 1/3 
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FQH-like state. We conjecture that the low-energy effec- 
tive theory for this fully spin-polarized ground state is 
characterized by the K matrix 



K 



fl 
-3 



with the filling fraction 

v = Q T K~ 



2/3. 



(2.9a) 



(2.9b) 



Clearly, this K matrix does not obey the decomposi- 



tion (1.2), since time-reversal symmetry is spontaneously 
broken. The degeneracy |deti^| = 3 is confirmed by the 
numerical results. The state thus obtained resembles the 
conventional double-layer 2/3 FQH state, with the dif- 
ference that the electron spins are not fully polarized. 

Case A = 1, U/V = 0: Possible paired state. A time- 
reversal symmetric state with a spectral gap and a three- 
fold ground state degeneracy is obtained for small U/V . 
This state cannot be captured by the Abelian Chern- 
Simons theory in Eq. since its degeneracy is not 

the square of an integer, despite the time-reversal sym- 
metry. One may speculate that this state realizes some 
real-space pairing of spin-up with spin-down electrons, 
since for small U/V it costs little energy to have two 
electrons of opposite spin at the same lattice site. 

We close Sec. |n] with some technical material. To de- 
termine the degeneracy of the ground state unambigu- 
ously, we have used spin-dependent twisted boundary 
conditions along the x direction defined by 



•0 


= 











(2.10) 



where \l/ is the many-body state, tr 3 acts on the spin de- 
grees of freedom, and x is the corresponding basis vector 
of the lattice. Imposing this boundary condition is equiv- 
alent to inserting the flux j x and its time-reversed flux 
— 7 X for electrons with spin up and spin down quantum 
numbers, respectively. This is a well-defined operation 
due to the residual spin-1/2 U(l) symmetry that pre- 
serves time-reversal symmetry for any value of j x . For 
the three cases discussed above, Figs.[l|b)-(d) show that 
the states of the (nearly degenerate) ground state mani- 
fold remain well separated from the continuum of states 
as 7 X is varied from to 2-7T, thereby confirming the nine- 
fold and three-fold degeneracy, respectively. 



III. EDGE THEORY WITH TIME-REVERSAL 
SYMMETRY 

We consider an interacting model for electrons in a 
two-dimensional cylindrical geometry as is depicted in 
Fig. [2j We demand that (i) charge conservation and 
time-reversal symmetry are the only intrinsic symmetries 
of the microscopic quantum Hamiltonian, (ii) neither are 
broken spontaneously by the many-body ground state, 




y 



+ Ly/2 



FIG. 2. Cylindrical geometry for a two-dimensional band 
insulator. The cylinder axis is labeled by the co-ordinate y. 
Periodic boundary conditions are imposed in the transverse 
direction labeled by the co-ordinate x. There is an edge at 
y — —Ly/2 and another one at y — +L y /2. Bulk states have 
a support on the shaded surface of the cylinder. Edge states 
are confined in the y direction to the vicinity of the edges 
y = ±L y /2. Topological band insulators have the property 
that there are edge states freely propagating in the x direction 
even in the presence of disorder with the mean free path I 
provided the limit t/L y <C 1 holds. 



and (iii), if periodic boundary conditions are assumed 
along the y co-ordinate in Fig. [2] then there is at most 
a finite number of degenerate many-body ground states 
and each many-body ground state is separated from its 
tower of many-body excited states by an energy gap. 
Had we relaxed the condition that time-reversal symme- 
try holds, the remaining assumptions would be realized 
for the FQHE. 

In the open geometry of Fig. [2j the only possible exci- 
tations with an energy smaller than the bulk gap in the 
closed geometry of a torus must be localized along the 
y co-ordinate in the vicinities of the edges at ±L /2. If 
L y is much larger than the characteristic linear extension 
into the bulk of edge states, the two edges decouple from 
each other. It is then meaningful to define a low-energy 
and long-wavelength quantum field theory for the edge 
states propagating along any one of the two boundaries 
in Fig. [2j which we take to be of length L each. 

The low-energy and long-wavelength effective quan- 
tum field theory for the edge that we are going to con- 
struct is inspired by the construction by Wen of the chi- 
ral Luttinger edge theory for the FQHEP^^l. As for the 
FQHE, this time-reversal symmetric boundary quantum 
field theory has a correspondence to the effective time- 
reversal symmetric bulk topological quantum-field theory 
built out of 2N A belia n Chern-Simon fields and defined 
by Eqs. (O) and (OM 



The simplest class of quantum Hamiltonians that ful- 
fills requirements (i)-(iii) can be represented in terms of 
2N real- valued chiral scalar quantum fields <P i (t, x) with 
i = 1, . . . , 2N that form the components of the quantum 
vector field &(t,x). After setting the electric charge e, 
the speed of light c, and H to unity, the Hamiltonian for 
the system is given by 



H :=H n 



(3.1a) 
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where 



:= I dx —dA r V dA , (3.1b) 

47T 



with V a 2N x 2N symmetric and positive definite ma- 
trix that accounts, in this bosonic representation, for the 
screened translation-invariant two-body interactions be- 
tween electrons. The theory is quantized according to 
the equal-time commutators 



64(4, a;') ^-■itt[k- ) 1 sgn{x-x') + Q 



where if is a 2N x 2N symmetric and invertible ma- 
trix with integer-valued matrix elements, and the ma- 
trix accounts for Klein factors that ensure that charged 
excitations in the theory (vertex operators) satisfy the 
proper commutation relations. We review the construc- 
tion of the vertex operators in detail in Appendix [A] In 
particular, fermionic or bosonic charged excitations are 
represented by the normal ordered vertex operators 



a T (x) = —a_ T (x). A discussion of the gauge symme- 
tries of this model and the properties of L can be found 
in Appendix [A] 

The theory ( |3.1[ ) is inherently encoding interactions: 
The terms H and H int encode single-particle as well as 
many-body interactions with matrix elements that pre- 
serve and break translation symmetry, respectively. Re- 
covering the single-particle kinetic energy of N Kramers 
degenerate pairs of electrons from Eq. (3.1b) corresponds 



to choosing the matrix V to be proportional to the unit 
2N x 2N matrix with the proportionality constant fixed 
by the condition that the scaling dimension of each elec- 
tron is 1 /2 at the bosonic free-field fixed point defined by 
(3.1c) Hamiltonian H . Of course, to implement the fermionic 



statistics for all 2N fermions, one must also demand that 
all diagonal entries of K are odd integers in some basis 
(see Ref. 125]) . While we shall proceed with the inter- 
acting bosonic theory here, we complement this analysis 
with a review of the stability of TV Kramers pairs of non- 
interacting fermions presented in Appendix [B] 



Time-reversal symmetry of the edge theory 



(3.1d) 



where the integer-valued 2A r -dimensional vector T deter- 
mines the charge (and statistics) of the operator. The 
operator that measures the total charge density is 



Z7T 



(3.1e) 



where the integer-valued 27V-dimensional charge vector 
Q, together with the if-matrix, specify the universal 
properties of the edge theory. The charge q T of the ver- 
tex operator in Eq. (3. Id) follows from its commutation 



with the charge density operator in Eq. (3.1e), yielding 
q T = T T Q. 

Tunneling of electronic charge among the different edge 
branches is accounted for by 



H iBt := - / da; ^ h T (x) : cos (t t K $(2) + a T (xfj : . 

TeL 

(3.1f) 

The real functions h T (x) > and < a T {x) < 2w en- 
code information about the disorder along the edge when 
position dependent. The set 



L := {T e 1 2N \T T Q = 0} , 



(3-lg) 



The operation of time-reversal on the <t fields is defined 

by 

T$(i,a;)r _1 := S x $(-*, x) + nR- 1 S ; Q, (3.2a) 
where 



- ( I J I and v 





J- ~ lo 1 



(3.2b) 



This definition ensures that fermionic and bosonic vertex 



operators as in Eq. ( 3. Id ) are properly transformed under 
time reversal. More precisely, one can then const ruct a 



pair of fermionic operators and of the form (13. ldh 



by suitably choosing a pair of vectors T x and T 2 , respec- 
tively, in such a way that the operation of time-reversal 



maps into +^2 whereas it maps \& 2 into Thus, 
it is meaningful to interpret the block structure displayed 
in Eq. ( 3.2b ) as arising from the upper or lower projection 



,t 



along some spin-1/2 quantization axis. 



Time- reversal symmetry on the chiral edge theory (3.1 ) 



demands the following conditions, which we explain be- 
low: 



V = +Y, 1 VY, 1 , 



(3.3a) 



encodes all the possible charge neutral tunneling pro- 
cesses (i.e., those that just rearrange charge among the 
branches). This charge neutrality condition implies that 
the operator &Jp(t, x) is bosonic, for it has even charge. 
Observe that set L forms a lattice. Consequently, if 
T belongs to L so does — T, In turn, relabeling T to 



-T in H int implies that h T (x) 



-h_ T (x) whereas 



K = -^K^, 



Q = 2iQ, 



(3.3b) 



(3.3c) 
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The first two conditions - Eqs. (3.3a) and (3.3b) - follow 



h T (x) — T (x), 



(3.3d) 



from the requirement that H be time-reversal invariant. 
In particular, the decomposition ( |1.2[ ) of K follows from 
Eq. (|3.3b|) and K = K T . The third condition - Eq. (|3.3cb 



- states that the charge density is invariant under time 
reversal. In pa rticul ar, the decomposition (1.2) of Q fol- 
,(x) = f-a SiT (i) + ttT t S 4 q) mod 2tt. (3.3e) lows from Eq. |3.3c| ). Finally, TH int T^ 



H;, 



requires 



h T (x) cos (r r K$(t,x) + a T (x) 



TeL 



h T {x) cos (t t K$(t,x) +a T (x) 

TeL 

h T (x) cos (- (Sj T) T .ftT $(-t, a;) + a T (x) - 7rT T Q 

reL 

X! cos (-T T i\T $(-t, a;) + a Ei T (a:) - T) T Q 

reL 

^ /i EiT (z) cos (T T ^ *(-t, i) - a Ei T (z) + tt^T)" 1 " E 4 Q 



(3.4) 



r 



leading to the last two relations - Eqs. ( 3.3d I and ( 3.3e ) 



as the conditions needed to match the two trigonometric 
expansions. 

Disorder parametrized by h T (x) = +h_ T {x) and 
a T {x) = —a_ T (x) and for which the matrix T obeys 



and 



E 1 T = -T, 



T T £ i Q is an odd integer, 



(3.5a) 



(3.5b) 



cannot satisfy the condition (3.3e) for time-reversal sym 



metry. Such disorder is thus prohibited to enter H int in 



Eq. (3. If I, for it would break explicitly time-reversal sym- 



metry otherwise. Moreover, we also prohibit any ground 
state that provides exp ^iT T K &(t, x)j with an expecta- 
tion value when T satisfies Eq. ( |3.5[ ), for it would break 
spontaneously time-reversal symmetry otherwise. 



B. Pinning the edge fields with disorder potentials 



Solving the interacting theory (3.1 ) is beyond the scope 
of this paper. What can be done, however, is to identify 



those fixed points of the interacting theory (3.1 ) that are 



pertinent to the question of whether or not some edge 
modes remain extended along the edge in the limit of 
strong disorder h T (x) — )> oo for all tunneling matrices 
TeL entering the interaction (3. If) . 



This question is related to the one posed and answered 
by Haldane in Ref. l26lfor Abelian FQH states and which, 
in the context of this paper, would be as follows: Given 
an interaction potential caused by weak disorder on the 



edges as defined by Hamiltonian (3. If I, what are the tun- 
neling vectors TgL that can, in principle, describe rel- 
evant perturbations that will cause the system to flow to 
a strong coupling fixed point characterized by h T — > oo 
away from the fixed point H Q 1 (See Ref. |45]for an an- 
swer to this weak-coupling question in the context of the 
IQSHE and Z 2 topological band insulators.) By focusing 
on the strong coupling limit from the outset, we avoid the 
issue of following the renormalization group flow from 
weak to strong coupling. Evidently, this point of view 
presumes that the strong coupling fixed point is stable 
and that no intermediary fixed point prevents it from 
being reached. 

To identify the fixed points of the interacting the- 
ory (3.1) in the strong coupling limit (strong disorder 
limit) h T — > oo, we ignore the contribution H and re- 
strict the sum over the tunneling matrices in H int to a 
subset H of L (H C L) wi th a precise definition of 
that will follow in Eq. ( 
follows the strong-coup 



3.10). For any choice of H, there 



ing fixed point Hamiltonian 



:= 



dx^ h T {x) : cos (t t K$(x) 

n Tern 



(3.6) 



We conjecture that a fixed point Hamiltonian (3.6) is 
stable if and only if the set H is "maximal". The 
study of the renormalization group flows relating the 
weak, moderate (if any), and the strong fixed points in 
the infinite-dimensional parameter space spanned by the 
non- universal data V, h T (x), and a T {x) is again beyond 
the scope of this paper. 

The reader might wonder why we cannot simply choose 
H = L. This is a consequence of the chiral equal-time 

as emphasized by Haldane 



commutation relations (3.1c 
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in Ref. [33 that prevent the simultaneous locking of the 
phases of all the cosines through 

d x (t t K $(t, x) + a T (x)^ = C T (x) , (3.7) 

for some time independent real- valued function C T (x). 
Even in the strong-coupling limit, there are quantum 
fluctuations as a consequence of the chiral equal-time 
commutation relations (3.1c) that prevent minimizing 



the interaction H int by minimizing separately each con- 



tribution to the trigonometric expansion (3. If I. Finding 



the ground state in the strong coupling limit is a strongly 
frustrated problem of optimization. 

To construct a maximal set H, we demand that any 
T € H must satisfy the locking condition (3.7). Further- 



more, we require that the phases of the cosines entering 
the fixed point Hamiltonian (3.6) be constants of motion 



d x (T T K$(t,x)) ,H B 



. 



(3.8) 



To find the tunneling vectors T £ H, we thus need to 
consider the following commutator 



d x (T T K$(t,x)) , h T ,(x) cos(T' T K4>(t,x')+a T ,(x'] 



r 



-\2-kT t KT' h T ,{x) sin (T' T K^>(t, x') +a T ,(x') 



(3.9) 



and demand that it vanis hes. This is achieved if 
T T KT' = 0. Equation EM implies that any set 



II is composed of the charge neutral vectors satisfying 
T T KT' = 0. It is by choosing a set H to be "maximal" 
that we shall obtain the desired criterion for stability. 



C. Stability criterion for edge modes 

We presented and briefly discussed in the introduction 
(see Sec. [l| the criteria for at least one branch of edge 
excitations to remain delocalized even in the presence 
of strong disorder. Here we prove these criteria. The 
idea is to count the maximum possible number of edge 
modes that can be pinned (localized) along the edge by 
tunneling processes. The set of pinning processes must 
satisfy 



T T Q = and T T KT' = 0, 



(3.10) 



which defines a set H introduced in Sec. IIIB| (Note, 
however, that EI is not uniquely determined from this 
condition.) It is very useful to also define the real exten- 
sion V of a set H, by allowing the tunneling vectors T 



that satisfy Eq. (3.10) to take real values instead of inte- 
ger values. Notice that V is a vector space over the real 
numbers. We shall also demand that H forms a lattice 
that is as dense as the lattice L by imposing 



vnL 



(3.11) 



For any vector T G V, consider the vector KT. It 
follows from Eq. plpf that KT _L T',VT' G V. So if 
maps the space V into an orthogonal space V . Since K 
is invertible, we have V 1 - — K V as well as V = K^ 1 
and thus dim V = dim V . Since dim V + dim V < 
27V, it follows that dimV < N. Therefore (as could be 
anticipated physically) the maximum number of Kramers 
pairs of edge modes that can be pinned is N; if that 
happens, the edge has no gapless delocalized mode. 



Next, let us look at the conditions for which the maxi- 
mum dimension N is achieved. If dimV = dimV -1 " = N, 
it follows that V V 1 - = R 2N , exhausting the space of 
available vectors, and t hus in this case the charge vector 
Qe? 1 because of Eq. ( |3.10| ). Consequently K~ X Q G V, 
and we can construct an integer vector T \\ K~ l Q by 
scaling K~ X Q by the minimum integer r that accom- 
plishes this (which is possible because K~ Y is a matrix 
with rational entries and Q is a vector of integers). Such 
a vector T G H is written as 



+(«• 



A) 
A) 



Q 



(3.12) 



The existence of (k — A) 1 follows from det K ^ and 



detK= (~) N [det(« — A)] 2 . 
Using T, we construct the integer 



(3.13) 



i?:=-T T S i g, (3.14) 

which, as we argue below, determines if it is possible 
or not to localize all the modes with the N tunneling 
operators. Here we employ Eq. (3.3e), also noticing that 
Y>iT — —T, and write 



nR 



ttT' E 



I' 



— otj,(x) — a s 2*0*0 J mod 2n 

( — ctrp(x) — a_f,(x)J mod 2-7T 
mod 2ir , 



(3.15) 



where in the last line we used that ct T (x) = —a_ T {x) for 



-L all T e L. If T satisfies Eq. (3.15), then R must be an 



even integer. If Eq. (3.15) is violated (i.e., R is an odd 



integer) then T is not allowed to enter H int for it would 
otherwise break time-reversal symmetry [thus h^(x) = 
must always hold in this case to prevent T from entering 
H- mi ]- We therefore arrive at the condition that 
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• If the maximum number of edge modes are localized 
or gaped, then R must be even. 

A corollary is that 

• If R is odd, at least one edge branch is gapless and 
delocalized. 

It remains for us to prove that if R is even, then one 
can indeed reach the maximum dimension N for the space 
of pinning vectors. This is done by construction. Take 
all eigenvectors of Si with +1 eigenvalue. We can take 
(N — 1) of such vectors, all those orthogonal to Q; for 
the last one we take T. One can check that these N 
vectors satisfy Eq. ([3.10 1 with the help of S 1 A"S 1 = 
—K [listed in Eq. ( pUbf ] and of f || K~ X Q. Now, the 
(N — 1) vectors E X T = +T are of the form T T = {t T ,t T ), 
where we need to satisfy T T Q = 2t T p = 0. This leads 



to T T Si Q even, and then Eq. (|3.3e|) brings no further 



conditions whatsoever. So we can take all these (N — 
1) tunneling vectors. Finally, we take T as constructed 
above, which is a legitimate choice since R is assumed 
even and thus consistent with Eq. (3.15). Hence, we have 



constructed the N tunneling vectors that gap or localize 
all edge modes, and can state that 

• If R is even, then the maximum number of edge 
modes are localized or gaped. 

As a by-product, we see that it is always possible to 
localize along the boundary at least all but one Kramers 
degenerate pair of edge states via the (N — 1) tunneling 
vectors that satisfy Yj-^T = +T. Thus, either one or no 
Kramers degenerate pair of edge state remains delocal- 
ized along the boundary when translation invariance is 
strongly broken along the boundary. 



we probe the stability of the FQSHE described by H 
when perturbed by H int (ReP^). 

To answer this question we supplement the condition 
T T Q = on tunneling vectors that belong to L and H, by 
T T S = 0. By construction, S is orthogonal to Q. Hence, 
it remains true that EI is made of at most N linearly 
independent tunneling vectors. 

The strategy for establishing the condition for the 
strong coupling limit of H int to open a mobility gap for 
all the extended modes of H thus remains to construct 
the largest set EI out of as few tunneling vectors with 
T = — SjT as possible, since these tunneling vectors 
might spontaneously break time-reversal symmetry. 

As before, there are (N — 1) linearly independent tun- 
neling vectors wit h T = +E 1 T, while the tunneling ma- 
trix T from Eq. (3.12) must belong to any EI with N 



linearly independent tunneling vectors. 

At this stage, we need to distinguish the case 



from the case 



T 5 = 0, 



T T S ^ 0. 



(3.17a) 



(3.17b) 



In the former case, the spin neutrality condition (3.16b) 



holds for T and thus the stability criterion is unchanged 
for the F QSHE. In the latter case, the spin neutrality 
condition (3.16b) is violated so that H int is independent 



of any tunneling matrix proportional to T. Thus, when 
Eq. (3.17b) holds, as could be the case when k oc 1 n and 
A = say, the FQSHE carried by at least one Kramers 
pair of edge states of H is robust to the strong coupling 
limit of the time-reversal symmetric and residual spin- 
1/2 U(l) symmetric perturbation H int - 



D. The stability criterion for edge modes in the 
FQSHE 

What is the fate of the stability criterion when we im- 
pose the residual spin-1/2 U(l) symmetry in the model so 
as to describe an underlying microscopic model that sup- 
ports the FQSHE? The residual spin-1/2 U(l) symmetry 
is imposed on the interacting theory (3.1 ) by positing the 
existence of a spin vector S = — E x S G I? N associated 
to a conserved U(l) spin current. This spin vector is the 
counterpart to the charge vector Q = +5^ Q £ 1 2N . The 
condition 

S = -T, X S, (3.16a) 
is required for compatibility with time-reversal symmc- 



IV. 



SUMMARY 



try and is the counterpart to Eq. (3.3c). Compatibility 



with time-reversal symmetry of Q and S thus implies that 
they are orthogonal, Q T S = 0. If we restrict the inter- 
action (3. If) by demanding that the tunneling matrices 
obey 



T T S = 0, 



(3.16b) 



We have considered in this paper a subclass of time- 
reversal-symmetric fractional topological liquids without 
quantized charge and spin Hall conductance. These 
states can be viewed as "zero filling fraction" quantum 
Hall states, that are related to an Abelian Chern-Simons 
bulk theory. The states we considered depart from previ- 
ous constructions that place together two copies of FQHE 
systems, and as such they do not need to satisfy spin con- 
servation or display quantized spin Hall conductances. 

We have analyzed the stability of the edge theory as- 
sociated with this type of state, and obtained a discrimi- 
nant, the parity of an integer, that resolves whether there 
remains or not delocalized edge states in the presence of 
disorder. When the discriminant is even, there are no 
gapless edge modes. In contrast, gapless edge modes are 
protected by time-reversal symmetry when the discrim- 
inant is odd. These results contain as particular cases 
those that display a quantized FQSHE^H where the dis- 
criminant has a relation to the then quantized spin Hall 
conductance. 
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We have also presented a concrete lattice realization 
of a FQSHE. There have been numerous studies on the 
effect of strong interactions in time-reversal-symmetric 
systems with non-trivial topology at half filling^. In 
contrast, we have considered the effects of interactions 
at a partial filling of bands. At 2/3 filling of a lattice 
with 24 sites, exact diagonalization delivers a ground 
state with nine-fold degeneracy, which we interpret as 
a time-reversal-symmetric fractional topological liquid. 
We studied the stability of this phase by tuning parame- 
ters of the electron-electron interaction. In particular, we 
found a transition toward a phase of spontaneously bro- 
ken time-reversal symmetry, which is related to a FQHE 
at 1/3 filling. 

Let us remark that the lattice model presented in this 
paper is an example of fractionalization in two spatial 
dimensions without breaking time-reversal symmetry. It 
thus joins the ranks with the other known examples 
thereof constructed on lattices so far, that of electron- 
fractionalization in graphenc-like systems in Refs. |4"81 - 
I5T1 together with the lattice gauge theory presented in 
Ref. [531 that of the triangular lattice quantum dimer 
model in Ref. |S5J and that of the doubled chiral spin 
liquids in Ref. OjSJ 

We would like to close this paper by spotlighting some 
perspectives on the differences between the Z 2 topolog- 
ical band (weakly interacting) insulators and the time- 
reversal-symmetric fractional topological liquid states 
whose very existences are driven by interactions, in par- 
ticular as to the importance that one should associate 
with the bulk and boundary states. In the case of the 
non-interacting (weakly interacting) systems, the edge 
states play a disproportionally important role, in that 
the bulk states are just band insulators without any 
ground state degeneracy. On the other hand, the time- 
reversal-symmetric fractional topological liquids display 
quite rich bulk phenomena, including the possibility of 
fractionalized quasiparticles, regardless of whether gap- 
less edge modes survive or not. Fractionalized particles 
can be probed without looking at the edge: capacitive 
measurements in the bulM^J, for instance, have revealed 
fractionalized electrons in the bulk of v = 1/3 FQH 
states. From this perspective, fractional time-reversal- 
symmetric topological liquid states should, as one might 
expect, be much richer in content than Z 2 topological 
band (weakly interacting) insulators. 
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Appendix A: Chiral bosonic quantum theory 

In this Appendix, we review the construction of 2A^ 
Fermi-Bose and 2N quasi-particle vertex operators from 
the chiral bosonic quantum fields $ i (t, x), i = 1, . . . , 2N, 
that enter the time-reversal invariant quantum edge the- 
ory with broken translation invariance (3.1) and dis- 



cuss their universal properties. To this end, we con- 
sider the universal data (K, Q) entering the theory de- 
fined by Eq. (3.1) as opposed to the non-universal data 



[y,h T (x),a T (x) 

On the chiral bosonic quantum fields we impose the 
boundary conditions 



K i} x + L) = Kij $j(t, x) + 2mi i 



(Al) 



with G Z for all i = 1, . . . , 2A^. Together with the 
condition that the tunneling vectors T are restricted to 
have integer-valued components, this ensures that the 
Hamiltonian H is single- valued. 

The chiral nature of the bosonic quantum fields arises 
from demanding that the equal-time commutator 

d i (t,x) ) ^(t,a/)] = -mfK^sgnix-x^ + Qij) (A2) 



holds for any pair i, j — 1, . . . , 2N. Here, 



(A3) 



and the antisymmetric 2A^ x 2N matrix L is defined by 
(see Ref. [ 



hj = sgn(i - j) {K i3 + Q i Q j 



(A4) 



where sgn(O) = is understood. It then follows that the 
quadratic theory (3.1b) is endowed with chiral equations 



of motion. Finally, we need to impose the compatibility 
conditions 



(_)*« = (_)< 



1, 



,2N, 



(A5) 



in order to construct local excitations with well-defined 
statistics. 

Define for any i = 1, . . . , 2N the pair of normal-ordered 
vertex operators 



t 

q-p,i 



(t,x) 



and 



*Li(*.*) 



(t,x) 



(A6a) 



(A6b) 



respectively. For any i — 1, . . . , 2N, the quasiparticle ver- 
tex operator x) is multi valued under the trans- 



formation ( A13[) provided |det_ftT| > 1 in contrast to the 
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Fermi-Bose vertex operator ${_ bi (f,3;) which is always 



single valued under the transformation (A13) 



For any pair i,j = 1, • • • , N, the equal-time commuta- 
tor (3.1c) delivers the identities 



and 



(A7a) 



(A7b) 



(A8a) 



number operator [recall Eq. (All] 

L 

1 







(A9a) 



*,(*,i)-$,(t,0) 



with eigenvalue one, while the Fermi-Bose vertex opera- 
tor ^j_ b At, x) is an eigenstate of the conserved operator 



1 

2^ 



(A9b) 



with eigenvalue one for any i = 1, . . . , 2N. 



(A8b) 



respectively. Here, the quasiparticle vertex operator 
^q.p i(t, x) is an eigenstate of the conserved topological 



The permutation statistics obeyed by any pair i,j = 1, . . . ,2N of quasiparticle and Fermi-Bose operators follows 
from the application of the Baker-Campbell-Hausdorff formula, 



and 



(AlOa) 



(AlOb) 



when x x', respectively. We conclude that, for any x ^ a/, demanding that the 2N x 2 A matrix K and the 2A- 
component charge vector Q are integer- valued together with the compatibility condition ( A5 ) is required to obtain 
local excitations carrying the Fermi-Bose permutation statistics 



i 



(AlOc) 



Let us now deduce the connection between the charge 
vector Q, the conserved operators C i , and the total 
charge density operator Q that follows from integrating 
Eq. (3.1e) along the edge. The charge vector Q enters 



for i = l. 



H, 



explicitly the theory after coupling the 2A chiral scalar 
fields to an external vector gauge potential with the com- 
ponents A and A 1 through the minimal coupling. The 
minimal coupling consists in replacing the x derivative 
by the covariant derivative 



. , 2 A and adding the contribution 

L 

■ urrcnt : = / dx ^ A (V C 'x® 




(Allb) 



on the right-hand side of Eq. (3.1a). The theory is then 



(Alia) 



invariant under the pure U(l) electro-magnetic gauge 
transformation 

$^$ + Qx, Ai^A 1 -d aX - (Allc) 
We can now define the total charge operator by 

Q := Qi ^ . (A12) 
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It follows that the charge associated with the quasipar- 
ticle operator fJ_ Pji and with the Fermi-Bose operator 

^l_ h i is given by K~^Q^ and Q it respectively. 

By assumption, the integer-valued 2N x 2N matrix K 
is symmetric and invertible. Consequently, its inverse 
K~ x is also symmetric, but its matrix elements are ra- 
tional numbers whenever |detif| > 1. Observe that the 
model (3.1) is invariant under the transformation 



$(t,x) $(t,x) +2ttT„ 



(Al3a) 



for any T t £ M. 2N that is independent of space and time 
and such that 



(Al3b) 



for al l tun neling vectors Tel. The quantum Hamilto- 

nian (3.1) thus possesses an emergent global (U(l)) 
symmetry compared to the microscopic model. The 
set of all r ational-valued vectors T, that satisfy condi- 
tions (A13) is the lattice L* dual to the lattice L. When 



|det K\ > 1, the lattice L is a sublattice of the dual lat- 
tice L* that is generated by the quasiparticles carrying 
a unit topological charge. The ground state of Hamilto- 
nian (3.1) with the periodic boundary conditions corre- 



sponding to the geometry of a torus is then degenerate 
with the degeneracy |det-ftT| > 1, which is nothing but 
the volume of the unit cell of the lattice L in units of the 
unit cell of the dual lattice L* . 

Note that there is no unique way to define the dual 
pair L and L*. For example, we can use Eq. (A7a) to 
define the dual lattice L* to be b ased on the hypercubic 
lattice Z 2N and, in turn, use Eq. (Al3b) to construct L. 
Alternatively, we can use Eq. (A8 



to define the lattice 
to be based on the hypercubic lattice 1 2N and, in turn, 
use Eq. ( |A13b[ ) to construct the dual lattice L*. Either 
ways, the ratio between the unit cells of the lattices L 
and L* is |detif|. We have chosen the latter option. 



Appendix B: Non-interacting fermionic edge theory 
with time-reversal symmetry 

In this Appendix, we review a non-interacting 
fermionic theory to analyze the role of disorder on the 
edge of a two-dimensional band insulator when imposing 
time-reversal symmetry. We consider a cylindrical geom- 
etry as is depicted in Fig. [5J If the length L y of the cylin- 
der is much larger than the characteristic linear extension 
into the bulk of edge states, the two edges at y = ±L y /2 
decouple. The low-energy and long-wave-length random 
single-particle Hamiltonian that describes any one of the 
two edges, each supporting N Kramers degenerate pairs 
of electrons, is then given by 



U{x) := iv F ct 3 <gs t N d x + W(x) = T~ x U{x) T, 

3 



(Bla) 



where the operation of time-reversal for spin- 1/2 elec- 
trons is represented by 



T 



ier 2 K 



(Bib) 



(K represents the operation of charge conjugation). Here, 
we have introduced the unit 2x2 matrix a n and the three 



Pauli matrices a 1 , a 2 , and ct 3 . In view of Eq. ( Blb[ ), a 1 , 



<7 2) and <7 3 are to be interpreted as the generators of the 
spin- 1/2 algebra of the electrons. The matrix 1^ is the 
unit N x N matrix. The matrix elements of the N x N 
Hermitean matrices W„(x) with fi — 0, 1, 2, 3 must obey 



W (x) 
W 2 {x) 



-WS(x), 

-w;(x), 



W 1 (x) 
W 3 (x) 



-Wl(x), 



(Blc) 



for time-reversal symmetry to hold. Hence, they can be 
taken as random numbers obeying the white-noise and 
Gaussian distribution of mean 



((W lt )..{x))=v i 5 ij S i 



(Bid) 



and co-variance 



W, 



x s ^ s i x ~ x ') > 



(Ble) 

with i,j = l,...,N the flavor index that label the 
Kramers degenerate pairs of electrons. The length scale 
i is the mean free path within the Born approximation. 
The channel fi = 3 represents forward scattering and 
the mean v { for i = 1, . . . , N results, through a gauge 
transformation, in a flavor-dependent shift of the Fermi 
velocity v F . 

It was known from the studies of quasi-one-dimensional 
wires in the 1980's that the random potential W(x) local- 
izes all TV Kramers degenerate pairs of electrons when N 
is everP^. It was only realized with the seminal work of 
Ando and Suzuura in Ref. [5S] on carbon nanotubes that 
the case of an odd number N is (i) of physical relevance 
and (ii) only localizes (N—l) Kramers degenerate pairs of 
electrons, leaving one pair delocalized along the edge^. 
It then took two groundbreaking papers from Kane and 
Mele on graphene with spin-orbit coupling to make the 
deep connection that this absence of Anderson localiza- 
tion is t he es sence of a two-dimensional Z 2 topological 
insulatorESEl. 

Observe that the integer quantum spin Hall effect 
(IQSH E) c an be deduced from the following special case 
of Eq. (Bl |. If we demand that the symmetry condition 



a 3 n(x)a 3 = H(x) 



(B2a) 



also holds in addition to the time-reversal symmetry in 
Eq. (Bl). Condition (B2a) is then nothing but the resid- 



ual spin-1/2 U(l) symmetry generated by rotations about 
the spin-1/2 quantization axis cr 3 . Imposing the residual 
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spin-1/2 U(l) symmetry (B2a) on the disorder potential 



W(x) amounts to the restrictions 
W 1 (x) = W 2 (x) 



0. 



(B2b) 



for all x along the edge. Condition (B2b) removes all 



backward scattering channels from the single-particle dis- 
order potential W(x). The single-particle Hamiltonian 
rl(x) thus decomposes into the direct sum of two Hamil- 
tonians, each of which realizes an integer quantum Hall 
edge, but with opposite quantized Hall conductivity of 
magnitude N in units of e 2 /h. The difference between 
these quantized Hall conductivities is proportional to N 
and yields the quantized spin Hall conductivity 2N in 
units of e/(47r). 

In Sec. |III| we consider an interacting effective quan- 
tum field theory including local multi-particle interac- 



tions that break translation invariance. The only allowed 
underlying microscopic symmetries of this interacting ef- 
fective quantum field theory are charge conservation and 
time-reversal symmetry. The time-reversal symmetry in 
Sec. Ill in particular Eqs. (3.3d) and (3.3e|, are inher- 



ited from the following properties of the single-particle 
disorder potential in Eq. (Bla). If A is any complex- 



valued matrix, denote with abs (A) and arg (A) the ma- 
trices with the matrix elements given by the absolute 
values and phases of the entries in A, respectively. One 
then verifies that 



abs (E x W(x) EJ = abs W(z), 

arg (E x W(x) Ei) = na x <g> E - argW(x), 



(B3) 



where E is the N x N matrix with one for all matrix 
elements. 
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